o 






ffi 



X 



H 



VISCORESISTIVE MHD CONFIGURATIONS OF PLASMA IN 

ACCRETION DISKS 

GIOVANNI MONTANI^^^^t ^nd RICCARDO BENINI^** 

^ENEA - C.R. FRASCATI (DEPARTMENT F.P.N.) VIA ENRICO FERMI, 45 (00044), FRASCATI (ROME), ITALY 

^ICRANET - C.C. PESCARA, P. DELLA REPUBBLICA, 10 (65100), PESCARA, ITALY 

^DEPARTMENT OF PHYSICS (G9) "SAPIENZA" UNIVERSITA DI ROMA, 

PIAZZALE A. MORO, 5 (00185), ROME, ITALY 

*ICRA - INTERNATIONAL CENTRE FOR RELATIVISTIC ASTROPHYSICS 

tMONTANI@ICRA.IT 

*RICCARDO.BENINI(aiCRA.IT 



rrj Abstract. We present a discussion of two-dimensional magneto-hydrodynamics (MHD) 



configurations, concerning the equilibria of accretion disks of a strongly magnetized as- 
trophysical object. We set up a viscoresistive scenario which generalizes previous two- 

,S^ dimensional analyses by reconciling the ideal MHD coupling of the vertical and the radial 

t^ equilibria within the disk with the standard mechanism of the angular momentum trans- 

^ port, relying on dissipative properties of the plasma configuration. 

^ The linear features of the considered model are analjrtically developed and the non- 

"^ linear configuration problem is addressed, by fixing the entire disk profile at the same 

Ci order of approximation. Indeed, the azimuthal and electron force balance equations are 

' ' no longer automatically satisfied when poloidal currents and matter fiuxes are included 

^^ in the problem. These additional components of the equilibrium configuration induce a 

K^ different morphology of the magnetic fiux surface, with respect to the ideal and simply 

^^ rotating disk. 

(3 Accretion disks, Plasma physics, 2-d MHD PACS:97.10.Gz; 95.30.Qd; 52.30.Cv 

m 

[^ 1. Preliminaries 

O 

f^ Understanding the mechanism of accretion that compact objects show in the presence of 



lower dense companions, is a long standing problem in astrophysics ( Bisnovatyi-Kogan and 



Lovelace, 2001 ). In fact, while in absence of a significant magnetic field of the massive body 



the disk configuration is properly described by the fluido dynamical approach (Shakura 



1973 ), the situation becomes a bit puzzling when we deal with a strongly magnetized source. 
Cd Since over the last three decades, increasing interest raised about the mechanism of the 

angular momentum transport across the disk profile. The solution of this problem stands 
as settled down in the form of a standard model, at least in the case of small (intrinsic) 
magnetic fields. This context prevents the formation of high magnetic back-reactions inside 
the plasma disk and the gas-like approximation is predictive. We remark how it is known 



from the observations (Verbunt, 1982| ) that the accretion profile takes place often with 
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the morphology of a thin disk configuration. In this limit of a thin gaseous disk, the 
hydrodynamical equilibria underlying the accretion process are cast into a one-dimensional 



paradigm, (see among the first analyses of the problem (Pringle and Rees, 1972 Shakura 



1973: Lynden-Bell and Pringle, 1974)). In such a fluidodynamical scenario, the accretion 



mechanism relies on the angular momentum transfer as allowed by the shear viscosity 
properties of the disk material. The differential angular rotation of the radial layers is 
associated with a non-zero viscosity coefficient, accounting for diffusion and turbulence 
phenomena. Indeed, microscopic estimations of the viscosity parameter indicate that the 
friction of different disk layers would be unable to maintain a sufficiently high accretion 
rate, and therefore non-linear turbulent features of the equilibrium are inferred. 

However, when the magnetic field B of the central object is strong enough, the electro- 



magnetic back-reaction of the disk plasma becomes relevant (Ruffini and Wilson, 1975) 



As shown in Coppi (2005), the Lorentz force induces a coupling between the radial and 



the vertical equilibria, which deeply alters the local morphology of the system. In particu- 
lar, the radial dependence of the disk profile acquires an oscillating character, modulating 
the background structure too. The existence of such a coupling breaks down the one- 
dimensional nature of the problem and suggests a revision of the original point of view of 
the standard model. 
Furthermore, 



m 



Coppi and Rousseau (2006) it is discussed how the plasmas, charac- 



terised by a /3-parameter (i.e. the ratio of the thermostatic pressure p to the magnetic 
one B'^/Stt) close to unity, exhibit an oscillating-like mass density and the disk is decom- 



posed in a ring profile. Thus, the analyses in Coppi (2005); Coppi and Rousseau (2006) 



demonstrate that the details of the disk equilibria are relevant in establishing a crystalline 
local structure inside the disk. This ideal MHD result constitutes an opposite point of 
view with respect to the idea of a diffusive magnetic field within the disk, as discussed in 



Bisnovatyi-Kogan and Lovelace (2001). The striking interest in the details of this local 



disk morphology relies on the idea that jets of matter and radiation are emitted by virtue 
of the strong magnetic field and of the axial symmetry of the accretion profiles (see for 
instance 



Lynden-Bell (1996)). 



Such an ideal two-dimensional MHD analysis is pursued by neglecting the accretion rate 



of the disk in the leading order (see also the global analysis presented in Ogilvie (1997)) 



and thus without facing the question about the azimuthal balance of the forces acting on 
the electrons. In fact, in these works (Ogilvie, 1997: Coppi, 2005 Coppi and Rousseau| 
2006), toroidal currents (and toroidal matter fluxes) are addressed only. We analyze in 



some detail the implications that accounting for the azimuthal and electron force balance 
equilibrium equations has on the local profile of the disk. We pursue the setting up of these 



two additional equations in the same approximation limit of the analysis in Coppi ( 2005 ) 



and in Coppi and Rousseau (2006) and we arrive to establish a direct relation between 



the radial infall velocity and the azimuthal Lorentz force. This result relies on addressing 
a link between the turbulent viscosity coefficient and the resistivity one. In particular, 
in the limiting case /? — )• oo, we are able to fix the vertical dependence of the plasma 



configuration, differently from the issue in Coppi (2005), for the same limiting case. The 
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radial velocity we derive this way has an oscillating behavior in its radial dependence and 
therefore has a zero net radial average. This matter infall profile is able to provide a local 



non-zero accretion rate and therefore it reconciles the crystal structure, outlined in Coppi 



(2005), with a real accretion feature active in the disk. 



The paper is organized as follows. In subsection |1.1| we describe the basic equations of 



the stationary MHD. In subsection 1.2 the standard hydrodynamic theory of accretion disk 
is presented. In Section [2] the bidimensional accretion disk is discussed in the framework of 
MHD. Section [3] is devoted to an analysis of the azimuthal equation and on the balance of 
the force acting on the electrons along the tangential direction. The linear approximation 
of this bidimensional model is discussed in Section [4j while the non linear configuration is 
developed in Section [5J Concluding Remarks follow 

1.1. MHD of the Steady State. The behavior of a plasma embedded into a magnetic 
and a gravitational fields can be described by a MHD approach for a wide range of con- 
figurations. In the case of hydrogen- like ions and when a steady state holds, the system 
configuration is governed by the following set of equations (as written in Cartesian coordi- 
nates, and we note that repeated indices are intended as summed). 



(la) di (ev,) = 

(lb) evidiVi = -dip + di 



diX+ -^ijkJjBk 



(ic) »,8,r +^rSm = o, 

where the mass density e and the total (electron plus ion) pressure p are related via the 
temperature T, according to the equation of state p = 2KTe/m, where K denotes the 
Boltzmann constant and m is the proton mass. Here T> denotes the shear viscosity co- 
efficient, Ji the density current vector, Vi the plasma velocity field, while, x stands for 
the Newton potential due to the central body of mass Mg (the self-gravity of the plasma 
being negligible), i.e. x{x^) = GM-^ j ^J bijX^ xJ . Instead, the plasma acquires an internal 
electromagnetic structure whose back-reaction can be relevant. 

The behavior of the total electric field Ei and magnetic field Bi is described by the 
Maxwell equations standing in Gaussian units as 

47r 
(2a) eijkdjBk = — Ji 

c 

(2b) diEi = Airp 

(2c) eijkdjEk = => Ei = -di^ 

(2d) diBi = {) => Bi = eijkdjAk , 

where Ai denotes the potential vector, p the electric charge density, ^ the electric potential 
and the electric and magnetic fields are related by the MHD condition 

(3) El = eimnVmBn ■ 
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Recalling the relation 

(4) (-ijkf-ilm = ^jl^km — SjmSkl 



and by means of Eq.(2a) the Lorentz force rewrites as 



(5) ^eijkJjBk = ]- {-BidiBi + BidiBi) . 

c 47r 

The scheme we traced above completely fixes the steady plasma configuration, once as- 
signed the central object morphology. 

1.2. The Standard Model of Accretion disks. The characterization of an accretion 
disk is obtained by addressing an axisymmetric MHD configuration (described by cylin- 
drical coordinates r, <j), z) for the plasma surrounding the central object. However, some 
fundamental features of the accretion process are fixed by the fluidodynamical approach 



describing the matter infall through a gas disk profile, as outlined in Shakura (1973). We 
aim to fix how the viscoresistive approach is not well-grounded on a microscopical level and 
a new perspective is required, as based on typical features observed in laboratory plasmas 



( Coppi 2005 ; Coppi and Rousseau , 2006 ) . 

In the standard model of accretion, the configuration of an axisymmetric thin disk is 
determined by the fluidodynamical equilibria, which take place within the central gravita- 
tional fleld. 

The thin character of the disk depth allows us to simplify the configuration problem, by 
integrating over the vertical profile, so fixing an effective one-dimensional hydrodynamical 



problem. When averaged out of the vertical direction ( Bisnovatyi-Kogan and Lovelace 



2001), the radial equilibrium reduces to the condition that the angular velocity ui of the 



disk takes the Keplerian profile 



(6) a;(r) = uk = ]j -^:^ ■ 

This statement is equivalent to neglect the role played in the equilibrium by the radial 
pressure gradient, retaining the centripetal force exerted by the central object as the dom- 
inant effect. This Keplerian nature of the disk is well-grounded and significant deviations 
from such a behavior are expected in advective dominated regimes only. 

The vertical equilibrium corresponds to the gravothermal configuration confining the 
disk profile, and it is therefore governed by the equation 

(7) —+u}^Z€ = 0, 

which, for an isothermal disk of temperature T and equation of state p = v^e (where 
Vg = 2K-qT /m is the sound velocity), gives the exponential decay of the mass density over 
the equatorial plane value eo(^)! i-e. 

(8) D{z') = ^ = exp (~ 
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Here, H = \/2Vg/L0^ provides an estimation for the real half-depth of the disk, which, in 
principle would have an infinite vertical extension. 

The azimuthal equilibrium describes the angular momentum transport across the disk, 
by virtue of a viscous stress tensor component Trij, which enters the relation (first integral 
of the azimuthal equation) 

(9) Md{L-Ld) = -27rr^Tr^ , Tr^ = Vr^^. 

Here D must be thought as a turbulent viscosity coefficient, while L is the angular mo- 
mentum per unit mass. Moreover, Lj is a fixed value of the specific angular momentum 
and 

(10) Md = -27rr / {evr)dz = -2-KraVr , 

J-H 

is the mass accretion rate, associated to the radial velocity Vr <Q and to the surface mass 
density a = J_^ edz. Finally the continuity equation 

(11) 1 d{revr) ^ djev,) ^ ^ 

r dr dz 

once integrated over the vertical direction provides with the fundamental relation Mj = 
const > 0. 

Recalling that the model is dominated by the Keplerian feature uj ~ wk and observing 
that diOK/dr = —3uj}<i/2r, we arrive to the following form for the angular momentum 
transport versus the plasma turbulent viscosity 

(12) Md{L - Ld) = SvrPwKr^ . 

Plasma viscosity. This picture of the accretion mechanism is however affected by a dis- 
crepancy between theory and observations of real systems. In fact, the viscosity coefficient 
D in the disk is too small if estimated by the microscopic plasma (or atomic) structure. 
Its microscopic features are unable to account for the accretion rates observed in some 
astrophysical systems, like X-ray binaries. The observed accretion rates, as estimated by 
the increasing disk luminosity £d ~ GM^M^/ro, (rg being the radius of the central object). 



would require a larger value of P, (see the discussion in Shakura (1973)). The solution to 
this discrepancy is inferred in the turbulent behavior arising in the disk when the Reynold 
number is sufficiently large. Such turbulence in the disk plasma would then be responsible 
for the appearance of the large value of the viscosity coefficient. 
Since, by definition L = wr^, we can infer 

(13) V = 2(TVtH/2, , 

vt being a turbulence velocity, given by vt = avs, where a is a free parameter. It must 
be noted that the axisymmetric disk is linearly stable with respect to small perturbations 
preserving its symmetry and the angular momentum conservation. For a discussion of the 



onset of turbulence by MHD instabilities, see the analysis presented in Balbus and Hawley 
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(1998). Such a review work is based on the Velikhov approach of 1959 and it requires 



non-hnear interaction among perturbations of very small amplitude. 



1.3. Comparison with the literature. The difficulty of this turbulent scenario, led B. 



Coppi to provide in Coppi (2005) a local 2-dimensional MHD formulation of the disk profile, 



which argues how the notion of plasma turbulence could be replaced by microscopic ring- 
like structures. In this framework, that we address below, the effect of turbulent viscosity 
has to be replaced by fundamental plasma instabilities, observed in laboratory experiments, 
like the so-called Resistive Ballooning Modes. An interesting ideal two-dimensional MHD 



analysis of a rotating disk equilibrium has been provided by G. I. Ogilvie in Ogilvie (1997), 



where the thin disk configuration is addressed in both the weak and strong magnetic field 
limit (also solutions are derived assuming self-similarity in the radial profile of the non- 
thin case) . The hypotheses at the ground of this work are the same as in Coppi ( 2005 ) ; 



m 



Coppi and Rousseau (2006), but the analyzed regimes are significantly different. In fact, 
Ogilvie (|1997) the radial dependence of the disk profile has a power- law structure and 



the existence of the ring morphology does not emerge for a strongly magnetized thin disk. 
The main reason of such significant deviation between these two ideal MHD approaches 
relies both on the local character of the analysis of Coppi (compared to the global profile 
of Ogilvie) and on the nature of the asymptotic expansion of Ogilvie, which privileges the 
second vertical derivatives with respect to the radial ones (so altering the Laplacian term 
appearing in the local equilibria fixed by Coppi). 

Our analysis differs from both these two ideal formulations of the axisymmetric MHD 
equilibria, not only for including dissipative effects (like viscosity and resistivity), but, 
overall in view of considering poloidal currents and matter fiuxes. In fact in order to 
describe the accretion phenomenon it is necessary to deal with radially in-falling material 
and expectedly a non-zero azimuthal Lorentz force. These two features are not present 



in the approaches followed in Ogilvie (1997); Coppi (2005); Coppi and Rousseau (2006) 



and our generalization involves dissipative effects too in order to compare the vertical and 

radial coupling, within the disk configuration, to the paradigm of the standard approaches. 

Indeed, the standard model treatment of the accretion process is properly reviewed by 



Bisnovatyi-Kogan and Lovelace (2001), where the main features of the disk morphology 



and the main successful issues are summarized in some detail. For a more recent analysis of 
the questions concerning the diffusive magnetic field living in the plasma configuration and 
its possible enhancement to get jet configuration see Spruit (2008). The models presented 



in these two reviews differ deeply from the analyses by Coppi or by Ogilvie, because the 
validity of the viscoresistive MHD approach is postulated on a phenomenological ground, 
almost disregarding the coupling between the radial and the vertical equilibria. 

The main merit of our approach is in reconciling the two-dimensional MHD scenario, 
including microstructures of the plasma configuration, with the the viscoresistive formu- 
lation of the standard theory. Despite we are here not able to solve the resistivity puzzle 



(see the dedicated paragraph of section 3.1), we demonstrate that the microstructures de- 
termined by B. Coppi for the ideal MHD approach still survive in the presence of non-zero 
viscosity and resistivity coefficients and of the related radial and vertical matter velocity 
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components. Demonstrating such a structural stability of the crystalline profile, we open 
the perspective to unify the microscopical and observational points of view, expectedly 
on the level of plasma instabilities. Anyway, the non-zero local accretion rate predicted 



by our model represents a significant improvement of the analysis in Coppi ( 2005 ) ; Coppi 



and Rousseau (2006), because it upgrades the crystalline picture toward a realistic disk 



configuration. 

2. TWO-DIMENSIONAL MHD MODEL FOR AN ACCRETION DiSK 

Let us now specify the steady MHD theory, discussed in Section 2, to the specific case of 
an accretion disk configuration around the compact (few units of Solar mass) and strongly 
magnetized (a dipole-like field of about 10^'^Gauss), i.e. a typical pulsar source. The 
gravitational potential of the pulsar has the form 

(14) x{r,z) = ^^, 

It is worth noting that the axial symmetry prevents any dependence on the azimuthal angle 
(j) of all the quantities involved in the problem. In this respect the continuity equation takes 



the explicit form (11), which provides the following matter flux associated with the disk 
morphology 

(15) ev = — dzQcr + eu)ir , z )res H — drQcz , 

r r 

where 0(r , z) is an odd function of z in order to deal with a non-zero accretion rate, i.e. 

(16) Md = -27rr / evrdz = 4-K@{r ,H)>0, 

J-H 

The magnetic field, characterizing the central object, takes the form 

(17) B = — dziJer + -es + -drijez, 

with ip = Tp{r , z^) and I = I{'ip , z). The similarity of the magnetic field and matter flux 
structure, is due to their common divergence-less nature. 

We now develop a local model of the equilibrium, as settled down around a radius value 
r = ro, in order to analytically investigate the effects induced on the disk profile by the 
electromagnetic reaction of the plasma. To this end we split the mass density and the 
pressure contributions as e = e(ro, z^) + e and p = p(ro, z'^) + p, respectively. The same 
way, we express the magnetic surface function in the form t/j = ^o('^o) + "01 ('^o > '^ ~ ^o > -2^), 
with Tpi <^ ipQ. The quantities e, p and ipi describe the change induced by the currents 
which rise within the disk imbedded into the external magnetic field of the central object. 
In general these corrections are small in amplitude but with a very short scale of variation. 
Thus, we are lead to address the "drift ordering" for the behavior of the gradient amplitude, 
i.e. the first order gradients of the perturbations are of zero-order, while the second order 
ones dominate. 
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As ensured by the corotation theorem Ferraro (1937), the angular frequency of the disk 



rotation has to be expressed via the flux function as OJ{^l>). As a consequence, in the present 
split scheme, we can take the decomposition w = wk + f^o^i, where li;k is the Keplerian 
term and u'q = const.. This form for uj holds locally, as far as (r — tq) remains a sufficiently 
small quantity, so that the dominant deviation from the Keplerian contribution is due to 

■01. 

Accordingly to the drift ordering, the profile of the toroidal currents rising in the disk, 
has the form 

c 



(18) 



J6 



AirrQ 



(a,Vi + 5,Vi) 



On the other hand, the azimuthal component of the Lorentz force is related to the existence 
of the function /(0, z) and it can be written as 

(19) ^0-7-^ {djdr^ - drld,^) . 

47rrQ 

2.1. Vertical and Radial Equilibria . We now fix the equations governing the verti- 
cal and the radial equilibrium of the disk, distinguishing the fluid components from the 



electromagnetic back-reaction (as develop by Coppi (2005); Coppi and Rousseau (2006)). 
Such a splitting of the MHD equations for the vertical force balance gives 



(20) 



D{z') =. 



exp 



(21) 

Here, as in Section 



eo('^o) 

dzP + oj^ze 



eo(?'o) = e(?'o, 0) 



mojt 



1 



Airrl 



[dliJi + a,Vi) d,^i 



0. 



1.2 



the behavior of the function D{z'^) accounts for the pure ther- 
mostatic equilibrium standing in the disk when the vertical gravity (i.e. the Keplerian 
rotation) is sufficiently high to provide a confined thin configuration, while the tempera- 
ture T admits the representation 

(22) 2KbT = m- = m- 



p + p 



2KBif + T). 



e e + e 

The radial equation underlying the equilibrium of the rotating layers of the disk, can be 
decomposed into the dominant character of the Keplerian angular velocity plus an equation 
describing the behavior of the deviation 5uj, i. e. 

(23) 



CJ ~ OJK + (^^W ~ tJo(V'o) + '^oV'l I 

2a;Kro(e + e)a;oV'i + - — 2 (^^^i + '^rV'i) ^rV'i '■ 



(24) 



47rrg 



dr 



p + 



Svrrg 



(drlpl 



+ 



AttTq 



dri)idl'il)i 



Here we neglected the presence of the poloidal currents, associated with the azimuthal 
component of the magnetic field. 
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Let US define the dimensionless functions Y, D and P, in place oi ipi, e and p, i.e. 

(25) y.Mi, D.^, P./3^, 

OroVo eo Po 

wliere po = 2K-QTeQ and /3 = 8itpq/Bq^ = l/(3e^) = /cg-ff^/S. Here we introduced the 
fundamental wavenumber ko of tlie radial equilibrium, defined as ko = Suj'^/v'j^, with 
v\ = i?2o/4vreo, recalling that B^o = Bz{r,z = 0) = dr^ilJo/rQ. Thus, we introduce the 
dimensionless radial variable x = fco(r — rg), while we assume that the fundamental length 
in the vertical direction be A = ,Je^H, leading to define u = z/A. By these definitions, 
the vertical equilibrium can be restated as 

(26) 5„2P + e,D-2 {dlY + e^dlv) d^2Y = , 
while the radial configuration is fixed by 



(27) 



(d + ^d\ Y + {d^Y + e.dlY) (1 + 9,y) + ^d^P = . 



The two equations above provide a coupled system for P and Y once the quantities D 
and D are assigned; hence we are able to determine the disk configuration induced by 
the toroidal currents. However, it remains to be described the accretion features of of the 
disk and the angular momentum transport across the crystalline structure lacking in the 



original analysis developed in Coppi (2005); Coppi and Rousseau (2006). 



3. AZIMUTHAL EQUATION 

The equilibrium along the toroidal symmetry of the disk is described by an equation 
which takes the exact expression 

(28) eVrdr{ojr) + eVzdz{ojr) + eujVr = -^dr {Vr^dr^) + dz \Vdz{'jor)] + F,^ , 

Fff, being the c/i-component of the Lorentz force. The corotation theorem, stating uj = cj{'ip), 
allows us to rewrite the equation above as 

ervrdrtp + ervzdztp + 2evr— = 



-dr [Vr^Jdr^) + —dz [Vruo'dz 






{oj' = dcj/dip). By virtue of the magnetic field form, we restate the l.h.s. of this equation 
in the form 

er^ {VrBz — VzBr) + 2eVr—f = 
(30) ^ 

-^dr {Vr^u:'dr^) + -^dz [Dru'dz^] + ^ . 
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In the local model we are addressing, by virtue of Eq.(19) the azimuthal equation stands 
at vq as 

erl {VrBz - v^Br) + 2eVr—j- = roPo (^rV'l + ^zV'l) + 

LUr, 

(31) 1 

where Pq = ^(''o)- This equation accounts for the angular momentum transport across 
the disk, by relaying on the presence of a viscous feature of the differential rotation in the 
spirit of Section 3. Furthermore, in this scheme, we link the radial and vertical velocity 
fields with the poloidal currents present in the configuration. Indeed we aim to get an 
equilibrium picture in which the disk accretion is induced by the poloidal currents directly, 
i.e. we want to establish a relation between Vr and F^. To this end, as well as for the 
model consistence, we have to analyze the structure of the electron force balance equation. 

3.1. The electron force balance equation. In the presence of a non-zero resistivity 
coefficient rj, the equation accounting for the electron force balance, reads as 

(32) E+- AB = r]J. 

c 

Since the contribution of the resistive term is expected to be relatively small, it turns out 
as appreciable only in the azimuthal component of the equation above. Thus, the balance 
of the Lorentz force has dominant radial and vertical components, providing the electric 
field in the form predicted by the corotation theorem, i.e. 

(33) E = — A S = — —Vip = (driller + S^V^) ■ 

c dip c 

Since the axial symmetry requires E^j) = 0, the (/>-component of the equation above takes 
the form 

(34) VzBr - VrB^ = C7]J^ . 

In the local formulation around rg, the azimuthal component stands as follows 

2 

(35) VrB, - v,Br = ^ {df.iPi + a>i) . 



We now observe that, substituting Eq.(35) in Eq.(|31|), we get 

(36) ^^ V ^ 4^ y 

This relation together with the electron force balance one ( |35[ ) provide a system for the 
two unknowns Q and /, when -01 (i.e. Y) and e (i.e. D and D) are given by the vertical 
(26) and radial (27) equilibria. 
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The resistivity puzzle. Before proceeding with our analysis, it is worth noting that the 
explanation for a high value of the resistivity coefficient of the disk plasma is necessary 
to make consistent the accretion scenario of the standard model, but it is hard to be 
provided and we are led to speak of a real resistivity puzzle. In fact, in the limit of 
a small electromagnetic back-reaction of the plasma, we deal with a negligible vertical 
dependence of the flux surface function (i.e. Br ~ 0) and the azimuthal component of 



the electron force balance Eq.(32) reduces to the simple form VrB^ = —rjJtj). Since in 
the considered limit, the toroidal current can not be significantly intense, the only way to 
ensure a sufficiently high radial velocity (able to account relevant accretion rates), consists 
of accounting on the role played by resistivity in fixing the equilibrium. However, like the 
viscosity coefficient, the quantity rj results to be very small if microscopically estimated and 
the question arises about the origin of strong resistive features of the disk. In the standard 
model, an "anomalous" resistivity coefficient is postulated in view of the turbulent behavior 
associated to the Velikhov instability of the rotating configuration. This scenario privileges 
a link between the viscosity and resistivity parameters, leading to values of order unity for 
magnetic Prandtl number. In the study we pursue below, despite the equilibrium is fixed 
for generic values of these dissipative coefficients, the linear and non linear solutions are 
derived under the assumptions of important viscoresistive effects. Our aim here is not 
to solve the puzzle of resistivity, but reconciling the crystalline profile with the presence 
of poloidal currents and matter fiuxes. For a proposal of the mechanism able to restore 



accretion in agreement to very small values of the resistivity coefficient, see Benini and 



Montani| ( |2009 ). 

4. The linear approximation 

Let us first analyze the linear model, corresponding to the request SroV'o ^ drijji and 
dzijji — 0. These conditions are equivalent to impose on the dimensionless vertical (|26|) 



and radial (27) equations, the restriction y <C 1 and to approach the limit e^ < 1 (i.e. 
j3 > 1). Furthermore the pressure gradient is neglected because it is not expected to be 
the responsible for a significant deviation from the Keplerian disk. Such an approximation 
will hold in the non- linear and low (3 values limits too. 

In the linear regime, we clearly have to require e = e-|-e~e = ^^{rQjD^z'^). Thus, the 
radial equation, in its linear form, reads as 

(37) dl^i + 5,Vi = -klD{z^)iji . 

In the same approximation, the linear electron force balance equation can be easily recast. 
In fact, neglecting the function D{z'^), as allowed by the much greater value of H with 
respect to A, at the lowest order, we get: 

(38) vr = -— _9.e = /^V'l , 

roe ATTdroWo 

We observe that (as in the dipole magnetic profile) drotpo < 0. 



In order to deal with the linear equation ( 36 ) which directly links the radial velocity Vr 
to the azimuthal Lorentz force Fj,, we are naturally lead to set the conditions rj = r]Q/D{z'^) 
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and Vq = (c^?7oeo)/(4vr) (r/o being the resistivity at ro), obtaining 



(39) 



SvrrgwKeo 



dzldroipo . 



Comparing the two expressions for Vr (p38|) and p9|), we get the compatibihty relation 
(40) 



2r]oeoc^{korofujK , 



(droi'Oy 



Equation ( 40 ) provides a relation between the two functions / and ipi . By ( 38 ) we easily 
get 

T/oeoc^fcoro 



(41) 



e 



47r 



dzY. 



Hence, we can find out the vertical and radial velocity components as 



(42) 



1 „ „ mc'^ko 



ro^o 



-drQ 



An 



dz{drY) , 



(43) 



-^ae = ^^5^y. 



^oeo 



47r 



where we neglected the function D[z^) because its contribution is here of higher order. 

Finally, we establish some phenomenological relations to shed light on the physical 
implications of the proportionality between the parameters Pq ^-nd tjq. According to the 
analysis developed in Section |1.2[ the turbulent viscosity coefficient can be expressed as 
follows 



(44) 



Po = -aeov^QH , 



where Vsq denotes the sound velocity on the equatorial plane. A reliable estimation for the 
resistivity is provided by the relation 



(45) 



% 



rUelAe 



where me is the electron mass, Ue the electron number density and z/ie a typical collision 
frequency between ions and electrons. Comparing equations (44) and (45 ) with the relation 
Vq = c^?7oeo/(47r), we can obtain the following estimation for H 



(46) 



H 



3 c^rUpUu 



Sir aneC^Vso Svr aUercVso 



where re is the classical radius of the electron. Equation ( 46 ) establishes a direct relation 
between the half depth of the disk and the collision frequency between electron and ions 
once the features of the background quantities Ue and Vso are assigned. 



viscoresistive mhd configurations of plasma in accretion disks 
5. Non-linear Configuration 
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In the non-linear regime, we have seen that the radial and vertical configurations (26 27 ) 
are described by the system in the unknowns Y and P, in correspondence to assigned D 
and D. In such a non-linear case, we retain also the direct relation between the turbulent 
viscosity coefficient and the resistivity one, which now reads 

AttVo SttoVsoH 



(47) 



V 



c2(e + e) c2(Z) + /3-iL)) ' 



reducing Eq. ( 36 ) to the searched link connecting Vr and F^ 

1 



(48) 



2uKeo{D + l3'Wy 



Substituting Vr by this equation into the electron force balance Eq.(35), we get the following 
expression for Vz 



(49) 



1 



dzMD + (3~W) 



n 



^ (9,Vi + dli^i] 



1 



{droi^a + dri^l) F^ 



eo ' ■ ' 2a;Keo 

We recall that, when addressing the splitting ip = ipQ + ipi, the azimuthal Lorentz force 
reads 

(50) F^ ~ 2 [9zl (SroV'O + dri^l) - drld^i^l] ■ 



In terms of the function 0, the expressions above for Vy. (48) and v^ (49) stand as 



(51) 

(52) 



dzQ 



2UJK 



dr& 



dzi^i 



Vo (5,Vi + S'V'i 



1 



2ujk 



{drolpO + dripl) F<t> 



The solution of the first equation can be taken as 



(53) 



G 



2WK 



dz{F^) , 



which, substituted in the second, yields an integro-differential relation for i<0, i.e. 
(54) 9,Vi / dzidrF^) = -2ujk'Do (9,^1 + 5,Vi) + (SroV'o + 5,Vi) F^ . 

Once the behavior of V'l is provided, i.e. Y{x, u), from the equation above we can determine 
the form of -F0(r, z) and eventually calculate the (/^-component I{ip, z) of the magnetic field. 
It is immediate to check that the integro-differential equation above provides, in the linear 
regime, the right compatibility condition we previously fixed in this limit, i.e. eq. (40). 



(55) 



Equation (54) can be easily rewritten in the dimensionless form 

duY I du{d^A^) = - [dlY + ezdlY) + (1 + d^Y) A^ 



A^ = F^/2uj-K^Voko being a dimensionless function. 
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5.1. Solution in the limit e^ — )• 0. The system constituted by the vertical (26) and the 



radial (27) equations admits a simple solution in the limit of vanishing e^, as outlined in 



Coppi (2005). In fact, such an asymptotic regime can be consistently represented by the 



following expressions for the magnetic surface Y and the pressure term P respectively 



(56) 



Y = F{u^)smx , P = Y^ = F^{u^) sin^ 



where F{v?) is a generic function of the vertical coordinate. In Coppi (2005) this function 
has been fixed by a supplementary condition, related with the higher order in e^. Now we 
show that in the present analysis F{v?') is naturally fixed by the azimuthal equilibrium, in 



particular by the partial integro-differential Eq.(55). It is immediate to check that if we 

set A(f,{u'^, x) = A{v? 

(57) 
(58) 



sinx, such an equation can be split into the two relations 

A^ = -Y 



dF 
du 



duF = F^ 



The solution of the ordinary integro-differential equation ( 58 ) can be taken in the form 

A exp < - 



(59) 



F u' 



} 



where A denotes a positive constant quantity. Thus, by Eq. (48), we get the following 
profile of the radial matter infall 



(60) 



-A exp 



{-\/^} 



smx 



2aVsokoH 
A exp 



-^A^^ 



smx , 



where we made use of the Shakura expression (44) for the viscosity coefficient Vq. The 



behaviour of the radial component of the velocity is depicted in fig{T] 

We see how the azimuthal equation completely fixes the disk profile at the zeroth order 
in e^. The resulting accretion radial velocity has an oscillatory character along the radial 



coordinate around vq, of the same kind found in Montani and Benini (2009). According 



with the discussion pursued in that work, this velocity profile is able to account for the 
local matter infall, but it does not clarify how the central object can receive this accretion 
rate. In fact the following expression for the accretion rate Mg holds 



Ms 



87r 



leading to the non zero average for x £ [— vr, vr 



Aa^/e^{x + koro) sinx , 



(62) 



/ 



Ms 
\ eoH'^Vso 



Wtt' 



Aa 



This behaviour is sketched in fig(2] 

The present analysis, differently from the linear case addressed in Montani and Benini 
(2009), lives in a full non-linear regime y ~ 1 (i.e. ^ ~ 1). Therefore we are able to 



identify the mass perturbation e, as the relevant feature to deal with a global matter infall. 
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Figure 1. We plot the radial component of the velocity Vr normalized 
to the constant sound velocity Wso in correspondence to the value a = 
OM,koH = 5,A= 1. 



Indeed, it seems that the non-linearity of the equilibrium configuration is not a strong 
enough feature to excite modes with non-zero radial averaged infall velocity. 

Finally we observe how the vertical decay of the magnetic flux surface, over the equatorial 



plane, is slower in the present case with respect to the analysis pursued in Coppi (2005), 
as outlined in Fi^ 

By other words, accounting for the azimuthal and electron force balance equations allows 
a higher intensity of the magnetic field in a given vertical region. The form of the total 
magnetic surfaces corresponds to the expression 

(63) 



Yt{x,u ) = X -\- Y{x,u ) = X -\- ^exp ( —Vu"^ ) sin(x) . 



Neglecting the toroidal component, which is regarded as small, the magnetic field compo- 
nents take the form 

/e^5„lT = ^/e^A^== exp (-Vv?) sin(x) , 



(64a) 



B 



20 
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Figure 2. The behaviour of the accretion rate as a function of the radial 
coordinate is depicted. It is outhned how the material in-falling from outside 
in tq is greater than the lost one. The values of the parameters for the dot- 
ted, dot-dashed, dashed and thin curve are given by a = 0.02, 0.02, 0.01, 0.03 
and e^ = 0.005,0.01,0.001,0.01 respectively. 
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Figure 3. The vertical dependence of the disk profile F(m^) is plotted in 
the present case (thin line) in comparison to the ideal MHD case (dashed 
line) . 



(64b) 



S, 



B 



zO 



dxX'Y = 1 + ^exp 



cos(x) , 
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The ^-field line are depicted in FigJ4J 
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Figure 4. The magnetic field lines neglecting the toroidal components, in 
correspondence to the values ^ = 1, ^/ej = 0.005 



6. Concluding Remarks 
Our analysis constitutes a bridge scheme between the standard one-dimensional model 



for an accretion disk ( Bisnovatyi-Kogan and Lovelace , 2001 ) and the reformulation provided 
in Coppi ( 2005 ) in terms of a local bidimensional MHD profile of the magnetically confined 



plasma. The original description relies on the setting of a viscoresistive MHD framework 
in which the angular momentum transport is allowed by the turbulent viscosity coefficient 
arising when the Reynolds number of the fluid has sufficiently high value. This paradigm 
has been able to provide with a satisfactory interpretation of some observed feature in 
strongly accreting system like X-ray binary pulsars. 



On the other hand, the approach introduced in Coppi (2005) starts from a criticism 



concerning the possibility to generate in the plasma of the disk significant viscosity and 
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resistivity effects. Tfierefore tfie disk morphology tias been addressed on a purely theoretical 
point of view, relying on the hypotheses of the ideal bidimensional MHD in axial symmetry. 
However this point of view cannot properly account for the accretion mechanism because 
it stands only for a velocity field completely dominated by the Keplerian rotation. In fact, 
including radial and vertical velocity in this bidimensional MHD scheme requires a proper 
setting of the electron force balance equation. 

We have constructed a scheme in which the viscoresistive bidimensional MHD scenario 



is implemented in the same spirit of Coppi (2005), but with the advantage to be able 



to account for the azimuthal and electron force balance equilibria, necessary to fix the 
accretion properties of the system. The resulting approach emerged as self consistent 
and, as far as a direct relation between the viscosity and the resistivity coefficients was 
addressed, we could fix the radial and vertical profile of the magnetic flux surfaces, due to 
the electromagnetic back-reaction. We devoted particular attention to the analysis of the 
linear case when the induced magnetic surface is a small correction only to the background 
features. Furthermore we solved the non- linear case in the limit of a vanishing ratio between 
the magnetic pressure to the thermostatic one. 

The results we derived above constitute the first step for a reformulation of the mecha- 
nism at the ground of the accretion phenomenon in terms of local structures arising in the 
disk morphology, instead of the diffusive description characterizing the standard model. 

This work was developed within the framework of the CGW Collaboration 
(www.cgwcollaboration.it). 
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